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Abstract 

Utilizing a recently derived extension of the Maxwell-Boltzmann 
distribution to low occupation numbers (Simons 1994) this investiga- 
tion discusses the structure of stellar dynamical isothermal spheres. 
The resulting models, which resemble King models, constitute a new 
sequence of physically well-motivated spherical equilibrium configu- 
rations of finite mass, and give nice fits to pre-core collapse globular 
cluster data. 



Subject headings: Galaxy: Globular Clusters: General; Celestial 
Mechanics, Stellar Dynamics; Galaxies: Elliptical and Lenticular, Cd; 
Galaxies: Fundamental Parameters. 



1 Introduction 



Globular clusters show lots of fine-structure related to mass-segregation, ve- 
locity anisotropies, core collapse, etc., but as a first approximation their 
structure (at least prior to core collapse) is surprisingly well described by 
simple King models based on truncated isothermal spheres (see Elson, Hut 
& Inagaki 1987 for a review). 

Globular clusters are small enough for two-body encounters among the 
stars to play a significant role in their evolution, so Maxwell-Boltzmann 
statistics (isothermal spheres) might be expected to be relevant for them, 
just like in the case of elastically scattering molecules in a gas. Elliptical 
galaxies are too large for two-body encounters among the stars to be im- 
portant, but if they formed through processes resembling violent relaxation, 
Maxwell-Boltzmann statistics could play a role there as well (see Binney & 
Tremaine 1987 for an overview and references on these topics). 

It has long been recognized, that a self-gravitating stellar dynamical sys- 
tem described by a Maxwell-Boltzmann distribution of the form f(E) oc 
exp(— (3E), where E is the stellar energy, has infinite mass, and therefore 
can at most provide an approximate description of, e.g., globular clusters or 
elliptical galaxies. 

In the context of model-building, this has been remedied in several ways 
in terms of introducing cutoffs in the distribution function. Most notably 
the King models (King 1966) based on the lowered isothermal sphere with 
fx(E) oc (exp(— f3E) — exp(— /3<E> )), where $o is a constant discussed below. 
King models fit the observational properties of globular clusters very well, 
and to less extent elliptical galaxies, but their somewhat ad hoc nature is 
slightly disturbing. 

Recently Simons (1994) derived an exact expression for Maxwell-Boltzmann 
statistics by a method which correctly accounts for discreteness effects by 
allowing macrocells in phase space to be occupied by integer numbers of par- 
ticles only. The resulting distribution, which forms the basis for the stellar 
dynamical applications in the present investigation, is given by 

F p =[Ag p exp(-f3E p )], (1) 

where [b] is the value of b rounded down to the nearest integer, F p and g p are 
the number of particles and the number of available states with energy E p , 
and A and (3 are constants to be determined from the conservation equations 
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J2 F p = N , mJ2 FpEp = E to t, where N and E tot denote total number of parti- 
cles and total energy, and m is the particle mass. In the present investigation 
all particles are assumed to be of equal mass, and the word "energy" is used 
as a short synonym for energy per unit mass, except E tot , which is a true 
energy. 

It can be shown explicitly (Simons, 1994), that this distribution maxi- 
mizes entropy under the implicit assumption of fixed energy levels. For self- 
gravitating systems no maximum entropy state exists, since the entropy can 
be increased by concentrating the core and letting some particles go into a dif- 
fuse halo. But in the case of e.g. globular clusters undergoing slow evolution 
dominated by two-particle encounters (i.e. prior to more violent phenomena 
like core collapse) the gravitational potential is almost time-independent, 
and stars can be expected to adjust into a distribution close to that which 
maximizes the entropy for a fixed potential. This is the distribution given in 
equation (1). 

In the following sections equation (1) will be discussed in the context 
of stellar dynamics, fits of globular cluster data will be obtained, and the 
fundamental parameters and assumptions in the model will be described. 

2 The Poisson equation for Simons- Maxwell- 
Boltzmann statistics 

To find the structure of a spherically symmetrical stellar system described 
by equation (1) one has to solve Poisson's equation on the form 



where the density is given by p = J f(E)d 3 v, and / denotes a mass per length 
cubed per velocity cubed; for stars of identical mass such a coarse-grained 
distribution is obtained from equation (1) by multiplying with i]/g p , where 
rj is the phase space density of a microcell occupied by one star. Macrocells 
will be constructed so as to contain the same number of microcells, g p . 

It is customary to introduce a relative gravitational potential, VP, related 
to the true potential $ via \P = — <3>+$o, and a relative energy e = — E+& = 
^ — \v 2 , with the constant $ chosen so that / > for e > 0, and / = 
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otherwise. Then ^ satisfies Poisson's equation on the form V 2x I/ = —4nGp 
with ^ — >• $o for r — > oo, or more explicitly 



with 



/(e) = 1 [exp(/3e)] . (4) 



This corresponds to <3> = \n{Ag p )/f3. 

1 /2 

Introducing a dimensionless radius s = r/a, where a = (j3 l l 2 g p /(Grj)} , 
and scaling potentials and energies in units of /3 one finally gets a dimension- 
less Poisson equation 



J = -16tt 2 jf [exp(/3e)](2(/3tf - /3e)) 1 / 2 d(/3e) 



s 2 ds \ ds 

jmax-l An(j+1) /-/3* 



-167T 2 £ / j(2(^-(3e)) 1 / 2 d((3e) + 3max (2(^ ~ /3e)) 1/2 d(/3e 

Jmax 

£ - lnj) 3 / 2 . (5) 



oll/2_2 Jmas 



Here j maa; = [exp(/3*)]. 

A solution of the dimensionless Poisson equation exists for each choice of 
scaled central potential, (3^(0), subject to the boundary condition (3^/'(0) = 
(no force at the center). The dimensionless potential decreases monotonously 
and approaches a negative constant, /3$o, for s — > oo. The corresponding 
density (in units of 77/ (g p (3 3 / 2 )) decreases from a central value, p*, and reaches 
when (3^/ = at the dimensionless tidal radius, s t . This is the major 
difference relative to the "traditional" isothermal sphere, where the density 
only asymptotically approaches zero. Thus the total mass, M = J 4npr 2 dr, 
is finite. 

To illustrate the importance of discreteness effects consider first the sim- 
plified situation where the physical gravitational potential as a function of 
radius, $(r), is fixed a priori (with $(r) ^ for r — > 00). In this case one can 
specify a partitioning of (x, v) phase space, including specification of energy 
levels, E p , with equal numbers of available states, g p . One may choose the 
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numbering of levels such that E\ is the lowest level, i.e. closest to $(0), etc. 
Energy level p is occupied by F p particles, where F p = [Ag p exp(— (3E P )] = 
[exp(/?e p )] = [exp(/?$ )exp(-/5E p )] = [exp(/3tf(0)) exp(/?($(0) - E p ))\, using 
the identities Ag p = exp(/3$ ), and e p = $ - E p = \1>(0) + $(0) - E p . 

The classical isothermal sphere has (fractional) occupation of energy lev- 
els all the way to E p — > for p — > oo (with particles reaching r — > oo), 
so to rescale the gravitational potential one would normally use $ = 0, 
\l/(r) = — $(r). For the true isothermal spheres of the present investiga- 
tion there is a finite number of particles, and levels are not fractionally 
occupied, so only levels with p below some finite p max contribute, corre- 
sponding to E p < E PmsiK < 0. Choosing $ = E PmSK the rescaled potential 
\P(r) = $ — $(r) now has the property that it is positive for exactly that 
range of radii where particles with $(0) < E p < E Pmix (^(O) > e p > 0) are 
able to move. 

The total mass M and physical energy (not per unit mass) E tot are given 

by 

M = mN = F p = m ^ [exp(/5^(0)) exp(/5($(0) - E p ))} (6) 

p p 

E tot = mY,E p F p = mJ2E p [exp(^(0)) exp(/5($(0) - E p ))\ . (7) 
p p 

If (3 were fixed, then a system with fixed total mass M could be realized 
for given particle mass, m, with that choice of ^(0) which obeys equation 
(0). Clearly, lowering m (and thereby increasing the number of particles, N) 
corresponds to picking a higher value of \]/ (0) to increase the sum. This is 
equivalent to choosing the (negative) value $ closer to 0, allowing particles 
to occupy states with higher "quantum number" p. Thus a decrease in parti- 
cle mass is directly related to an increase in the scaled potential, f3^, for fixed 
total mass. As one sees from the figures, solutions with higher dimension- 
less central potential (3^(0) (corresponding as just argued to more particles 
for fixed total mass, physical potential and phase space partitioning) are as 
expected closer to the standard isothermal sphere solution. I.e., discreteness 
effects are less conspicuous for large N. 

In reality, (3 would not be fixed in this procedure since the solution must 
also obey the equation for E tot , but since (for dimensional reasons) E tot = 
M (3~ x h(/3ty (0)) , where h is a slowly varying function of (3^(0), (3 does not 
change much if M and E tot are fixed. 
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For the self-consistent case where the potential $(r) is given by Poisson's 
equation rather than fixed a priori, the picture is more complicated, but 
$(r) is not strongly varying with N for fixed M and E tot , so qualitatively 
the argument above remains valid — discreteness effects are more pronounced 
for low N, and for fixed M and E tot (and partitioning of phase space) a 
change in m (N) corresponds to a change in "shape parameter" (3^(0) for 
the isothermal spheres. 



3 Comparing with observations 

As explained above, the new Maxwell-Boltzmann statistics in a natural way 
accounts for finite mass systems. However, many parameters appear (A, g p , 
f], (3), so it is not obvious that much has been gained in terms of constructing 
useful stellar dynamical models. It will now be demonstrated how these 
parameters are related to observable properties of, e.g., globular clusters. 

Figure 1 shows the density profiles for several choices of (3^(0). Densities 
are scaled in units of the central density, and radii in units of the King radius, 



r ° = V ^G P (0) ~ V 4^°- (8) 

Also shown in Figure 1 is the standard isothermal sphere, which the present 
models approach for (3^(0) — > oo, and the corresponding King models. The 
similarities are striking — not surprisingly, since King models also have an 
energy cutoff at E = $ and an element of rounding down (f K oc (exp(/3e) — 1) 
rather than [exp(/3e)]). 

Equally similar are the RMS velocities at a given spatial radius in Figure 
2. (For the full Simons-Maxwell-Boltzmann statistics one can show that 



The velocities shown in Figure 2 are yv 2 (r)(3/3 ). 

Figure 3 shows the so-called concentration parameter defined as c = 
log 10 (r</ro). The value of c is typically 0.3 units lower than the corresponding 
King model with the same (3^(0). King models fit globular cluster data quite 
well for 0.75 < c < 1.75 and elliptical galaxies moderately well for c > 2.2 
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(Binney & Tremaine 1987; Mihalas & Binney 1981). Similar limits are rea- 
sonable for the isothermal spheres. This corresponds to 5 < (3^(0) < 9 and 
P^(0) > 11 for the isothermal spheres. 

Each choice of dimensionless central potential, (3^(0), leads to values 
for the dimensionless central density, p* = p(0)g p (3 3 / 2 /r/, the dimension- 
less total mass, M* = MG/3/a, the dimensionless tidal radius s t = r t /a, 
and the dimensionless, central line-of-sight velocity dispersion, a* 2 = [3a 2 = 
2 J? Ei(*(a) " In j) 5/2 ds/5 - In jfl 2 ds. 

Explicitly expressing observable quantities in terms of i] (in units of 
M pc~ 3 (km/s)~ 3 ), (3 (in units of (km/s)~ 2 ), and g p , one gets 

P(0) = p^P^Me pc- 3 , (10) 

M = 3.546 x 10 3 M*(3- 3/4 r]' 1/2 gl /2 M , (11) 

r t = l^s^-^gy 2 pc, (12) 

and 

a 2 = a* 2 /?" 1 (km/s) 2 . (13) 

Notice that 7] and g p always occur in the combination i]/g p (g p is further 
coupled to A via <3> — l n (A^ p )//3, but this can be considered as a consistency 
equation for A). 

A typical model with (0) = 8, that gives a good fit to the surface 
brightness structure of globular clusters, has s t = 0.153, p* = 4.677 x 10 4 , 
M* = 0.1290, a* 2 = 0.9650. For M = 6 x 1O 5 M and a 2 = (7km/s) 2 this 
gives f3~ l = (7.126km/s) 2 , rj/g p = 2.10 x 10~ 4 M Q pc~ 3 (km/s) -3 , correspond- 
ing to a tidal radius of 60.3 pc and a central density of 3.6 x 10 3 M o pc~ 3 . 
Choosing rj to be higher than the maximal (central) phase space density of 
the system (rj > 10M Q pc~ 3 (km/s)~ 3 ) means that g p < 5 x 10 4 . Thus statis- 
tical mechanics gives reasonable results for a globular cluster of 6 x 10 5 1M 
stars if one divides phase space in microcells smaller than w 0.1pc 3 (km/s) 3 
and macrocells containing of order 10 4 or less microcells. 

4 Discussion and conclusions 

It has been demonstrated in this investigation, that a spherical stellar dy- 
namical equilibrium configuration described by Maxwell-Boltzmann statis- 
tics, contrary to the "standard" isothermal sphere, has finite mass. This is 
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because phase space cells can be occupied by an integral number of stars 
only, rather than the continuos distribution implicit in the classical Maxwell- 
Boltzmann distribution. The complete Simons-Maxwell-Boltzmann distribu- 
tion (Simons 1994) explicitly takes this into account, and in stellar dynamical 
applications this leads to finite mass because one does not integrate fraction- 
ally occupied cells over infinite space. 

A new class of physically well-motivated models for globular clusters were 
introduced based on the new distribution function. They resemble King 
models because of some common features such as an energy cutoff and a 
lowering of the classical Maxwellian, but in the King models these features 
are not included in a quite correct manner (King models still allow fractional 
occupation of phase space cells, and the lowered isothermal sphere does not 
maximize entropy, which equation (1) does for fixed gravitational potential, 
total energy and mass, resembling the situation for globular clusters prior to 
core collapse). 

The new models presented here fit the overall properties of globular clus- 
ters prior to core collapse well for approximately the same choice of central 
potential as for King models. They have constant "temperature" (/5 _1 ), but 
in contrast to the infinite mass isothermal sphere, this does not reflect itself 
in constant velocity dispersion (cf. Figure 2). For large central potentials the 
models get closer to the standard isothermal sphere, and the velocity disper- 
sion gets closer to (5. For very low central potential (j ma x = [exp(/3\I/)] = 1) 
the models go over into n = 3/2 polytropes, but this is far below the range 
of potentials giving reasonable fits to globular clusters. 

Any use of distribution functions in stellar dynamics involves an (often 
implicit) assumption about the "graininess" of phase space and a choice of 
unit cell. Contrary to quantum mechanics, where Planck's constant defines 
a fundamental cell size, there is no natural choice in stellar dynamics. In 
fact there are even several "schools" arguing either for a choice related to 
the observational resolution (interpreting entropy as a measure of informa- 
tion content), or to the typical phase space size of elements which are not 
disrupted by interactions over the lifetime of the system studied. The pa- 
rameters extracted from fits to observations for the present model include 
measures of the fundamental phase space density, 77, and the number of mi- 
crocells in a macrocell, g p , but only in the combination i]/g p (which can also 
be interpreted as the mass of a fundamental particle (star) divided by the 
volume of a macrocell). This could indicate, that there is a physically moti- 
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vated relation between the central potential (or more generally the structure 
of the cluster) and the scale on which mixing in phase space takes place. As 
discussed above, values giving fits to globular cluster data seem reasonable, 
and are not very dependent of the choice of model, as long as (3^(0) is not 
too small. 

I thank Jens Hjorth for useful discussions. This work was supported by 
the Theoretical Astrophysics Center under the Danish National Research 
Foundation. 
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Figures 



Figure 1. Density profiles of isothermal spheres with central potential 
(3^(0) = 3, 6, 9, 12 (solid curves), and oo (dotted curve; equivalent to 
the classical, infinite mass isothermal sphere). For comparison are shown the 
corresponding King models for the same choices of /3^(0) (dashed curves). 

Figure 2. One-dimensional root-mean-square velocity dispersion profiles for 



Figure 3. The concentration parameter c = log 10 (r t /r ) as a function of the 
dimensionless central potential for isothermal spheres (solid curve) and King 
models (dashed curve). 



the same models depicted in Figure 1. wrms = yv 2 (r){3/3. 
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